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For turbulence, although the two-point velocity difference u(x + r) — u[x) at each scale r has been studied 
in detail, the velocity average [u(x + r) + u(x)]/2 has not thus far. Theoretically or experimentally, we find 
interesting features of the velocity average. It satisfies an exact scale-by-scale energy budget equation. The 
flatness factor varies with the scale r in a universal manner. These features are not consistent with the 
existing assumption that the velocity average is independent of r and represents energy-containing large-scale 
motions alone. We accordingly propose that it represents motions over scales > r as long as the velocity 
difference represents motions at the scale r. 



I. INTRODUCTION 

Suppose that a velocity component u(x,t) is obtained 
at time t along a one-dimensional cut x of a turbulent 
flow. The ensemble average (it) has been subtracted, so 
as to have (u) = anywhere below. Then, the two-point 
velocity difference u_ and average u + are defined at each 
scale r as 

u—(r, x, t) = u(x + r, t) — u(x, t), (la) 

u(x + r,t)+u(x,t) 
u + {r,x,t) = - . (lb) 

While u- has been studied in detail^ u+ has not thus 
far. We are interested in it + and are to study its statis- 
tical features. 

Sreenivasan and Dhruva 3 studied the probability den- 
sity distribution of u+(r) in atmospheric turbulence. The 
shape of the distribution was close to that of it, regardless 
of the scale r. 

Tatsumi and Yoshimura 4,5 studied the probability den- 
sity distribution of it+(r) at each scale r in a theoretical 
manner, by applying a closure approximation to homo- 
geneous isotropic turbulence. 

Hosokawa 6 found an exact relation for turbulence that 
is homogeneous in the x direction. Since (u(x + r) 3 ) — 
(u(x) 3 ), where (•} denotes an ensemble average, Eq. (TTJ 
yields 

(u 3 _(r)) = -12{ul(r)u_(r)). (2) 

Through this equation, Kolmogorov's 7 - four-fifth law for 
the longitudinal velocity at small scale r in the inertial 
range 

(u 3 _(r)) = ~(e)r (3a) 

is equivalent to 

(u\{r)u-{r)) = ±-{e)r. (3b) 
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Here e is the energy dissipation rate. Kholmyansky and 
Tsinober 8 confirmed Eq. (|3b[) in atmospheric turbulence. 
Exact relations were also found by Hill 9 and Germano.^ 

These studies tend to assume that, although it_(r) rep- 
resents motions at the scale u+(r) does not depend 
on r and represents energy-containing large-scale motions 
&lone£&&& If this were the case, (it^(r)it_(r)) ^ for 
r in the inertial range [Eq. (|3bj) ] would imply a corre- 
lation between small-scale motions in the inertial range 
and large-scale motions in the energy-containing range. 6 ' 8 
Such a correlation is of interest because still controver- 
sial is whether or not small-scale motions are statistically 
independent of large-scale motions 

However, it is necessary to reconsider the assumption 
for it+(r). We find that u 2 , satisfies an exact scale-by- 
scale energy budget equation that involves Eq. (|3b|) . We 
also use experimental data to find that (u\(r)) / ' (u 2 + (r)) 2 
varies with r in a universal manner. Thus, albeit domi- 
nated by large-scale motions, u+(r) reflects motions at 
the scale r££. We accordingly propose an assumption 
that u+(r) represents motions over scales > r. Since the 
correlation of u 2 + (r) with it_(r) in Eq. (|3b[) is attributable 
to motions at the scale r, there is no more need to invoke 
a correlation between small- and large-scale motions. We 
discuss this and other implications. 



II. EXPERIMENTAL DATA 

Experimental data of grid turbulence, boundary layer, 
and jet are used here. They were obtained with a wind 
tunnel in our recent worki 1 ^ We outline the experiments 
and the data. Their parameters are shown in Table HI 



A. Grid turbulence 

For the grid turbulence, we set a grid across the wind 
tunnel. The grid was made of rods with 0.04 x 0.04 m 2 
in cross section and 0.20 m in separation. At a; w t = 3.5 m 
downstream of the grid, a hot-wire anemometer was used 
to measure the streamwise velocity U + u(t wt ) and the 
spanwise velocity v(£ w t)- Here U is the average, whereas 
w(iwt) and v(t wt ) are temporal fluctuations. The data 
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TABLE I. Parameters of grid turbulence (GT), boundary layer 
(BL), and jet (J) taken from Ref. 15: kinematic viscosity v, mean 
rate of energy dissipation (e) , Kolmogorov velocity ux , rms veloc- 
ity fluctuations (u 2 ) 1 ' 2 and {v 2 ) 1 / 2 , skewness factor (u 3 ) / (u 2 ) 3 ' 2 , 
flatness factors (it 4 ) / (u 2 ) 2 and (v 4 ) / {v 2 ) 2 , Kolmogorov length rj, 
correlation lengths L u and L v , and Reynolds number Re>, for the 
Taylor microscale A = [2(v 2 ) / {{d x v) 2 )} 1 / 2 . 



Quantity 




Units 


GT 


BL 


J 


V 




cm 2 s~ 1 


0.143 


0.143 


0.139 


(e) = 15v((d*vf)/2 




2 -3 

m s 


2.81 


12.6 


2.60 


u K = (^<e)) 1/4 




m s _ 1 


0.0796 


0.116 


0.0776 


(« 2 } 1/2 




m s~ 1 


0.696 


2.37 


1.56 


{v 2 ) 1 ' 2 




m s~ 1 


0.683 


1.96 


1.36 


<« 3 >/<« 2 > 3/2 
<« 4 >/(« 2 ) 2 






+0.08 


-0.10 


-0.04 






3.00 


2.69 


2.60 


(v 4 )/{v 2 ) 2 






2.98 


3.05 


3.05 


„ = (, 3 /<e» 1/4 




cm 


0.0180 


0.0123 


0.0179 


Lu — /q 00 (u(x + r )u(x)) dr / (u ) 


cm 


17.5 


43.0 


128. 


L " = f£°i v ( x + r ) v (* 


))dr/(v 2 ) 


ecu 


4.46 


5.68 


10.2 


Re A = (v 2 )^ 2 X/ V 






285 


1103 


1183 



are long so that their statistical significance is high. Also 
at x w t — 3.25m and 3.75m, we measured U + u(t w t). 
This measurement was not simultaneous with the above 
measurement but was under the same condition. 

The temporal fluctuations w(t w t) and v(t w t) measured 
at position x w t are converted into the spatial fluctuations 
u(x) and v(x) at time t, by using Taylor's frozen-eddy 
hypothesis 

x = -Ut wt and t = -jj 1 - (4) 

Since the grid turbulence was stationary in i wt and was 
decaying in x wt , it is homogeneous in x and is decaying 
in £ i 14 i 16 The data at x wt = 3.25 m and 3.75 m are used 
to estimate dt (it 2 } and so on at the time t corresponding 
to x w t = 3.5 m. 



B. Boundary layer and jet 

The data of the boundary layer and jet were obtained 
in the same manner as for the grid turbulence. For the 
boundary layer, roughness was set over the floor of the 
wind tunnel. The measurement was done in the log-law 
region at a streamwise position x wt where the bound- 
ary layer had been well developed. For the jet, a nozzle 
was set within the wind tunnel. The measurement was 
done where the flow had become turbulent. Since these 
measurements were under stationary conditions, we use 
Taylor's hypothesis to obtain spatial fluctuations that are 
homogeneous in the x direction [Eq. ([I])]. 

To explore common features, we compare the boundary 
layer and jet with the grid turbulence. The comparison 
is based on the correlation length L u (Table |T| , which 
represents large scales and also serves as the typical size 
of energy-containing eddies. 



III. MEANING OF VELOCITY AVERAGE 

Usually, statistics of u_(r) and v_(r) are discussed by 
assuming that u_(r) and v~{r) represent motions at the 
scale ri 2 - To be exact, we also have to consider motions 
at scales < r. Their contributions to u_(r) and u_(r) are 
nevertheless negligible as long as r is not too large ^ This 
is because the energy distribution of a turbulent flow is a 
continuously and rapidly increasing function of the scale 
r up to about the correlation length L u . 

Then, as long as r is not too large for (r) and w_ (r) 
to be safely assumed as representatives of motions at the 
scale r, we assume u+(r) and as representatives of 

motions over scales > r. If this is the case, u+(r) and 
v + {r) reflect motions at the scale r, albeit dominated by 
motions at the larger scales, > r. Those at scales < r also 
contribute to u + (r) and v + (r), but they are expected to 
be negligible as long as they are negligible in u_(r) and 
v_(r). 

Our assumption is justified by using the experimental 
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FIG. 1. Statistics for grid turbulence as a function of r/rj. (a) 
(u 2 + )/(u 2 ), {v 2 + )/(v 2 ), (u 2 _)/u 2 K , and (v 2 _)/u 2 K . The circles de- 
note {u 2 )/(u 2 ). (b) Flatness factors {u 4 + ) / (u\ ) 2 , (v\) / {v 2 + ) 2 , 
{ut_) / (ul_) 2 , and {v± )/{lP_ ) 2 ■ (c) Energy budget in Eq. {TO)}. The 
left-hand side ("tot.") and the second term in the right-hand side 
("trs.") are normalized by —3dt{u 2 )/2. The first term ("dis.") is 
normalized by 15v((dxu) 2 ). The dotted curve denotes the sum of 
the terms in the right-hand side. The arrows indicate r = L u and 
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FIG. 2. Statistics for boundary layer as a function of r/rj. (a) 
(u\)/{u 2 ), {v 2 + )/{v 2 ), {u 2 _)/u' 2 K , and {v 2 _)/u 2 K . The circles denote 
(u 2 )/(u 2 ). (b) Flatness factors (v\) / (v\) 2 and (v± }/{v 2 _ ) 2 . (c) 
Energy budget in Eq. mil l. The first term in the right-hand side 
("dis.") is normalized by 15v((d x u) 2 ), and the second term ("trs.") 
by 15v((d x v) 2 ) /2. The dotted curve denotes their sum. The arrows 
indicate r = L u and L v . 



data to compare (u+(r)} with (u 2 (r)). Here, 



u(r, x, t) = - 
r 



x-\-r 



u(X,t)dX, 



(5) 



which defines the motions over scales > r because those 
at scales < r have been smoothed away. The results are 
shown in Figs. QJa), Ufa), and|3](a). We observe {u 2 + (r)) ~ 
(u 2 (r)) for r < L u , where it is safe to assume tt+(r) as 
a representative of motions over scales > r. This is not 
the case for r > L u , where motions at scales < r are 
no longer negligible in (u^(r)) as well as in (u 2 _{r)) and 
accordingly these two are constant. 

Our assumption is self-consistent. For example, in the 
limit r — > 0, u + {r) and v+(r) reduce to u and v that 
represent motions over all the scales. 

The scales of motions could be also defined with wave 
number k. We have 

poo 

(u 2 _(r)) oc / [1 - cos(kr)]E u (k)dk, (6a) 
Jo 

poo 

{u%(r)) oc / [1 + coa(kr)]E u (k)dk, (6b) 
Jo 



where E u (k) is the energy spectrum. This is consistent 
with our assumption because 1 + cos(fcr) is maximal at 
k = and falls to zero at k = w/r, where 1 — cos(fcr) is 
maximal. 1? However, being localized in the x space, u+ 
or u_ is not localized in the k space. For the scales, our 
definition based on Eq. ([5]) is preferable. 

Phenomenologically, u- (r) and V- (r) represent the de- 
formation and rotation of an eddy with size r, respec- 
tively, if r does not exceed the typical size of energy- 
containing eddies, L u . This is just the r range for u_(r) 
and u_(r) to represent motions at the scale r [Figs.QJa), 
HI a )i Eta)]. Likewise, we assume u+(r) and v+(r) as the 
advection of an eddy with size r. The advection is due to 
sweeping by nearby eddies with sizes > r, being consis- 
tent with the above assumption that motions over scales 
> r are represented by u + (r) and v + (r). 

Between u + and u_, exact statistical relations exist 
if the turbulence is homogeneous in the x direction. 6 
We have already mentioned (v?_(r)) = — \2(u\(r)u- (r)) 
[Eq. ©]. In addition, since {u{x) 2 ) = {u{x + r) 2 ), Eq. {TJ 
yields 



(u 2 _(r)) 



(u+(r)u_ (r)) = 0. 



(7a) 
(7b) 



The same relations exist between v + and u_. For each 
r, the mean total energy (u 2 ) is decomposed into the 
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mean energies (u+(r)) and (v?_(r)) [Eq. (Tfal) ]. Although 
u+(r) does not correlate with tt_(r) [Eq. (17b)) ]. u^_(r) 
could correlate with u_(r) [Eq. {2}]. This correlation is 
through motions at the scale r, if we ignore motions at 
scales < r. 

The relations in Eqs. ([2]) and {7J) are not restricted to 
the homogeneous turbulence. The same relations exist in 
turbulence that is isotropic throughout the scalesi 8 ' 19 i 20 
We nevertheless focus on the homogeneous turbulence, 
which is more common than the isotropic turbulence. 
For example, an experiment is usually under a station- 
ary condition. The measured temporal fluctuations are 
converted by Taylor's hypothesis into spatial fluctuations 
that are homogeneous in the x direction [Eq. {2}]. 



IV. ENERGY BUDGET EQUATION 



By substituting Eqs. ([2]) and (ITal) into Eq. (|8a|) and also 
by rearranging the terms in the left-hand side, we obtain 



0_ 

di 



> 2 



15 



f (u 2 + {R))R 4 dR 
Jo 



30^ d{u\{r)) I5{u\{r)u-{r)) 
r dr r 



(8b) 



Thus, u 2 ^ satisfies an exact scale-by-scale equation. This 
would not exist if u+{r) were a representative of large- 
scale motions alone. It does represent motions at scales 
> r. 

To clarify the meanings of Eqs. (|8a|) and f|8b[) . their 
left-hand sides are rewritten by using the relation 



The scale-by-scale energy budget is studied in terms 
of u 2 + . We consider turbulence that is homogeneous and 
isotropic throughout the scales and also is decaying freely 
in time, for which an exact equation is known : 11 ! 21 



d_ 

In 



--(« ) H - 

2 V ' 4r 5 



f (u 2 _{R))R 4 dR 
Jo 

Vov d(u 2 _(r)) 5{u s _(r)) 



2r 



dr 



4r 



15 
15 



(u 2 _(R))R 4 dR 

f ' (u 2 + (R))R 4 dR 
Jo 



(9) 



which is from the asymptote (u(x + r)u(x)) — > for r 
oo. The results are 



d_ 

di 



15 

4r 5 



f (u 2 _(R))R 4 dR + f (u 2 _(R))R 4 dR 

Jo r->oo J 



l$vd{u 2 _{r)) _ 5(u 3 _(r)) 
2r dr 4r 



(10a) 



d_ 

dt 



15 

r 5 J,, 



XR))R 4 dR 



15 



.(i?))i? 4 di? 



30^ d(u 2 + (r)) 15(u%(r)u-(r)) 
r dr r 



(10b) 



These equations offer different descriptions of the same 
energy budget. By using Eqs. ([2]) and (|7al) . the left-hand 
side of Eq. (jlObl) and the individual terms in its right- 
hand side are able to be obtained one-by-one from those 
of Eq. (|10a|) . Thus, each two of these corresponding terms 
have exactly the same values. Since Eq. (|10a[) describes 
the budget of kinetic energy per unit time and mass in 
terms of u 2 _ (see below), it is natural to consider that 
Eq. (110b[) describes the same energy budget in terms of 
u 2 + . Here u 2 _ or u 2 + is not used as the energy itself but is 
used to describe its budget. Although u 2 _ and u+ behave 
differently [Figs. [Ua), [Ha), and[3Ja)], they are statisti- 
cally dependent on each other through Eqs. © and (Tfal) . 

For the grid turbulence, Fig. flTc) shows values of the 
terms of Eqs. (|10ap and (|10b[) . 22 Their corresponding 
terms are actually not distinguishable at all. 

The left-hand side of Eq. dTU|) is just a weighted average 
of d t (u 2 _(R)) or d t (u 2 y(R)) over R > r, which describes 
the mean rate of the total loss of the energy over scales 



> r. It reduces to — 39 t (u 2 }/2 = (e) as r — > and to as 
r — )■ cxd. The first term in the right-hand side of Eq. (p~0|) 
describes the mean rate of the energy dissipation into 
heat over scales > r. It reduces to 15^((<9;cu) 2 ) = (e) as 
r — > and to as r —> oo. These two are balanced with 
the second term, which describes the mean rate of the 
energy transfer across r from the larger to the smaller 
scales^ 

The transfer term for u 2 + involves u_ [Eq. (|10b[) ]. which 
is phenomenologically related to the typical time scale 
r/(u 2 _(r)) 1 / 2 for the energy transfer across the scale r 
(Sec. EJ. 

Between each two of the corresponding terms of 
Eqs. (|10a[) and (llObj) . there is a sign difference. This 
difference is associated with the fact that u 2 _(r) tends 
to increase with r while uz,(r) tends to decrease with r 
[Figs. [Ha), (Ha), and[3I»; see also Eq. {75}]. 

When the scale r is small enough to lie in the inertial 
or dissipative range, the left-hand side of Eq. (|10a[) is able 
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to be replaced by (e): 



I5v d(u 2 _(r)) 5(u 3 _(r)) 



2r 



Or 



4r 



(11a) 



This equation in turn applies to small scales of any tur- 
bulence that is isotropic over these small scales.-' 11 If the 
turbulence is also homogeneous in the x direction, Eq. 
(|10b() is likewise rewritten as 



<£> = - 



30vd(u 2 + (r)) 15(u 2 + (r)u_(r)) 
r dr r 



(lib) 



For the boundary layer and jet, Figs. (He) and [He) show 
values of the terms in the right-hand sides of Eqs. (|lla|) 

and $rn$M 

Within the inertial range, the right-hand sides of Eqs. 
(Illap and (jllbl) are dominated by the second terms. 
They correspond to Eqs. (|5aj) and (|3bl) . That is, while 
the transfer term for u 2 _ in Eq. (|lla[) corresponds to 
Kolmogorov's? four- fifth law (|5aj) , the transfer term for 



in Eq. (jllbl) corresponds to Hosokawa'sS. relation (j3b 
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FIG. 5. Same as in Fig. [4] but for boundary layer at r/rj 
10\ 10 2 , 10 3 , 10 4 , and 10 5 . 



V. ENERGY TRANSFER RATE 

To understand how the mean rates of energy trans- 
fer — 5(ui(r))/4r and 15(u+(r)u_(r))/r in Eqs. (JTDJl and 
are determined as observed in Figs. [TJc), [UJc), and 
[3fc), we study fluctuations of u_(r) and w+(r)u_(r). 

Figure Ufa) shows the probability density distribution 
of U-(r) in the grid turbulence. If r < L u ~ 1 x 10 3 ??, 
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FIG. 4. Statistics for grid turbulence at r/rj = 10 1 , 3 x 10 1 , 
10 2 , 3xl0 2 , 10 3 , 3x 10 3 , and 10 4 as a function of u~/{ui) 1/2 . 
(a) Probability density distribution of w_. (b) Conditional 
average (u 2 s r ) u _ normalized by The value at each r is 

shifted by 0.1. 
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the distribution is skewed with a long tail in the negative 
side. This tail explains the observed positive value of 
— 5(u^(r)}/4r. The value vanishes if r > L u , where the 
distribution of u_(r) is symmetric. 

Figure H{b) shows the average {u 2 + (r)) u _ conditioned 
on u-(r) in the grid turbulence. If r < L U) (u+(r)) u _ is 
greater at u_ = it, > than it is at w_ = — u* < 0. This 
asymmetry overbalances the negative skewness of w_ and 
explains the observed positive value of 15{vA_(r)u- (r)) /r. 
The value vanishes if r > L„, where (u+(r)) u _ is sym- 
metric. 
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The same r dependences are observed for the boundary 
layer in Fig. [5] with L u ~ 3 x 10 3 ?y and for the jet in Fig. [5] 
with L u ~ 7 x 10 3 r). 

Besides the asymmetry, (u 2 [ _(r)) u _ has features that 
are not common to the grid turbulence, boundary layer, 
and jet [Figs.|4|b),[5][b), andEJb)]. The configuration for 
turbulence production affects large-scale motions, which 
in turn affect vfr[r) at any r. However, not affected is the 
mean rate of energy transfer 15(ic_ (r)ii-(r))/r because 
correlates with u_(r) only through motions at the 
scale r (Sec. UlI)) . 

The local rate of energy transfer is known to fluctuate 
markedly This fluctuation could be related to those 
of — 5u 3 _(r)/4r and lhu\.{r)u-(r)/r, although they are 
not exactly identical to the local rate. 

Phenomenologically, when an eddy typical for the size 
r deforms and breaks apart into the smaller daughter ed- 
dies, the deformation energy (u^(r)) of the mother eddy 
is transferred to those of the daughter eddies. Since the 
time scale is r/(w?_(r)) 1 / 2 , this offers a conceptual expla- 
nation for the mean rate of energy transfer — 5(u 3 _(r))/4r. 
The mother eddy also has energy (tt5_(r)) for the advec- 
tion, which is also transferred to those of the daughter 
eddies. With the time scale r/(u?_(r)) 1 / 2 , this explains 
the mean rate in terms of 15(u?_ (r)u_(r))/r. These two 
explanations are equivalent if we consider eddies with all 
the sizes, since the advection u + of an eddy reflects the 
deformations it_ of the surrounding larger eddies. In fact, 
we have — 5(tt 3 _(r))/4r = lh(u 2 + (r)u^{r)) jr. 



These asymptotic values are the same if (u 4 ) / (u 2 ) 2 = 3, 
which is close to the case of the grid turbulence (Table 
HJ. However, in more detail, the flatness factors of u+(r) 
and v+(r) vary with the scale r, exhibiting maxima at 
r ~ L u and minima at r ~ 10~ 1 L U . Albeit of small 
amplitude, these variations are still above the statistical 
uncertainties that are seen as rapid fluctuations at r > 
10 1 L U . The same variation exists for the flatness factor 
of «+(?') in the boundary layer and jet [Figs. [2jb) and 

Hb)]^ 

For the observed variations of the flatness factors, we 
attempt an explanation, which is not complete but would 
serve as a guide for future studies. When r is not too 
large, u + (r) represents the advection of an eddy with 
size r, which is due to sweeping by nearby eddies. If each 
of them has the same sweeping direction, the correspond- 
ing ti_|_ value tends to lie at the tail of the u + distribu- 
tion. Since eddies with sizes < r do not contribute to 
the sweeping, an increase in r shortens the above tail. 
The flatness factor becomes small. With a further in- 
crease in r, eddies with sizes > r become sparse. The 
u+ distribution becomes enhanced at u+ ~ 0, i.e., cases 
of no eddies. The flatness factor becomes large. Finally, 
when r exceeds the typical size L u of energy-containing 
eddies, their random relative motions become dominant. 
The flatness factor tends toward the Gaussian value of 3 
[Eq. (HU)]. 



VII. CONCLUDING REMARKS 



VI. SCALE DEPENDENCE OF MOMENTS 

The r dependences of (u™(r)) and (u™(r)) are distinct 
from those of (u"(r)) and (i>"(r)). For example, while 
(u 2 _(r)) varies by several orders, {u 2 + (r)) varies only from 
(u 2 )/2 to (u 2 ) [Figs.EQ»,[2i;a), and^a)]. 5 This range is 
too narrow to exhibit a power-law scaling oc such as 
those known for (u™ (r)) 

The r dependences of (w" (r)) and (i>™ (r)) are also sub- 
ject to large-scale motions that depend on the configu- 
ration for turbulence production. For example, although 
exactly isotropic turbulence exhibits (u^i_(r)) — at all 
r, our grid turbulence exhibits some variation of (u+(r)). 
There still exist features that do not depend on the flow 
configuration, as demonstrated below by using long ex- 
perimental data. 

FigureQJb) shows the flatness factors in the grid turbu- 
lence, e.g., {iti(r)) / ' (u 2 + {r)) 2 . With a decrease in r, those 
of u_(r) and u_(r) increase as a result of the small-scale 
intermittency^ Those of u + (r) and u+(r) are constant 
within ±1%. The reason is the predominance of large- 
scale motions. In fact, 

{(u 4 ) 

, „.„ as r — > 0, 
{U) (12) 

2W + 2 ^ r ^°°' 



The velocity averages u + (r) and v + (r) represent mo- 
tions over scales > r, as long as r is not too large for 
motions at scales < r to be negligible, i.e., for the veloc- 
ity differences u_ (r) and v_ (r) to be safely assumed as 
representatives of motions at the scale r. Phenomenolog- 
ically, while U-(r) and u_(r) represent the deformation 
and rotation of an eddy with size r, and v + (r) 

represent the advection of such an eddy. 

The motions over scales > r are defined by u(r), i.e., 
velocity smoothed over the scale r [Eq. When r is 
not too large, we observe ~ (u 2 (r)) [Figs. Ufa), 

Ha), and Eta)]. 

We have found an exact scale-by-scale energy budget 
equation for u\ [Eq. (|10b|) or (jllbj) ]. which involves the 
relation found by Hosokawa 6 [Eq. (|3b[) ] as the mean rate 
of energy transfer across the scale r. This equation de- 
scribes the same energy budget as the equation for u_ 
[Eq. (jlOap or (| 1 1 a|) ] . Between them, the corresponding 
terms have exactly the same values. 

We have also found that the flatness factors of u+(r) 
and v+ (r) vary with r in a manner that is independent of 
the configuration for turbulence production [Figs. Htb), 
dtb), andEtb)]. Thus, albeit subject to large-scale mo- 
tions that depend on the flow configuration, u + and v+ 
have universal features. 

The analogues of u + and u_ exist in the orthonor- 
mal wavelet transformation, where u{x) is transformed 
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„27-29 



+00 



+00 +00 



(13a) 



m— — 00 



n=N m— — 00 

with two self-similar families of basis functions 



S„, ro (x) =2 n / 2 S(2 n x-m), (13b) 
W„ iTO (:e) = 2 n / 2 T^(2 n x - m). (13c) 

They are localized around x = 2 _ ™mr,, where r* is some 
unit scale. While S n , m (x) represents motions over scales 



> 2" 



Wn tm (x) represents motions at the scale 



r 

r ~ 2~™~ 1 rv If 2~ n r* is not too large, the transform 
uw is analogous to the velocity difference m_ . 27 ' 29 The 
transform us is analogous to the velocity average it+. 
There exist various wavelets based on various definitions 
of the scales. Those for Haar's wavelelj 2 ^— are the same 
as defined for the present study [Eq. (J5J)] . 

Tatsumi and Yoshimura 4,5 proposed a closure approxi- 
mation that u+ is statistically independent of it_ at each 
scale, which is interpreted as an approximation of statis- 
tical independence between the advection tt+(r) and de- 
formation u- (r) of eddies with each size r. This is exact 
in many cases [Eq. (I7b|) ] . An exception is the correlation 
of u 2 + with u_ studied here [Eq. ©J. 6 

There is a controversy as to whether or not small-scale 
motions in the dissipative and inertial ranges are statis- 
tically independent of large-scale motions in the energy- 
containing range £&&l±rM. To address this issue, u™ (r) at 
small r is often compared with u™(r) or u m , by assuming 
that they represent the large-scale motions. Being actu- 
ally a representative of motions over scales > r, u™(r) 
could correlate with it™ (r) through motions at the scale 
r [Eq. ©]. The situation is the same for u m , which corre- 
sponds to w™(r) for r — >• and hence could correlate with 
w"(r). Such correlations are spurious as those between 
small- and large-scale motions. In particular, the large- 
scale motions are not responsible to (u+(r)u~(r)) 7^ 
for r in the inertial range [Eq. (|3b[) ]. 

To conclude this paper, we underline that further stud- 
ies of u+ and v+ are promising^— ^r— For the above and 
other issues, u + and v + could reveal new features and of- 
fer new insights. This is because the description of turbu- 
lence by u + and v + differs from the existing description 
by u_ and V-, and also because u+ and v+ satisfy exact 
equations like Eq. (llObp that serve as the firm theoretical 
basis. 
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